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Abstract 


We give a quantitative theory of line forraatlon in an atmosphere 
above a surface with hackscattering properties. Sufficiently high spatial 
and spectral resolution spectra of resonance lines in lo region A can yield 
data on the surface scattering properties as vrell as the number density of 
scattering molecules. 

We discuss raacroscopically homogeneous models of scattering from the 
surface of lo and conclude that multiple reflection from crystal facets Is 
the most likely cause for the observed phase variations of the gewnetric 


albedo. 



1. Introductian. 

lo has a notably high albedo for visible and near Infrared radiation. 
Between 0.5 oatd 3 pa the geometric albedo of the leading side of :he satellite 
is observed to range from 0.7 to 1.1 (Johnson and McCord, 1970; Johnson, 1971; 
Johnson and Pilcher, 1974) in contrast to the value 0.67 for a laabert sphere 
with conservative scattering. 

Brown ^ (1975) have pointed out that, for such a high albedo, the 

sodium D-lines might be observed in absorption, if Region A* alone could be 
observed. 

This effect is Important as regards the search for other atomic resonance 
lines in the atmosphere of lo, and in the first part of this paper we present 
a quantitative theory. It is apparent from this work that observations of the 
sodium D-lines in region A could yield information about the photometric 
properties of the surface of lo as we3.1 as the density of sodium atenns in 
its atmosphere. 

In the second part of this paper we give a speculative but quantitative 
theory for the observed photometric properties of lo's surface. The high 
geometric albedo implies a scattering function with a stronger backward 
scattering component than for a lambert surface. This backward scattered 
component is seen in the phase variation of the albedo at 0,55 pm, which 
decreases by about 30% from 0® to 12* phase angle (Morrison et al. , 1974). 

In the first part of the paper we treat this component as a narrow bundle of 
constant intensity, centered on the incident beam. 


*Reglon A Includes the visible disk and the atmosphere gravitationally bound 
to the satellite. Sodium D-lines in emission have been observed in regions more 
distant from the satellite. No spectra from Region A alone are available at the 


present time. 
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Chanced backward scacter Is understood for dark surfaces such as lunar 
oaterlal (Hapke, 1963). As pointed out by Oetking (1966) these theories, 
which Involve shadowing effects, cannot be applied to high albedo surfaces, 
such as lo. Oetking showed that nufflerous white materials (e.g. , HgO) have 
enhanced backscatterlng; others, on the contrary, have enhanced scattering at 
the specular angle (Barkas, 1939). As far as we are able to judge, Ifeman 
(1955) first pointed to the possiblity of cube comers or "cat's eyes" asso- 
ciated with surface crystalline material. He was concerned with the observed 
negative polarization at small phase angles exhibited by many planetary sur- 
faces, but the same idea can obviously explain enhanced backscatterlng. Al- 
though the exact nature of the reflecting surfaces is unspecified in our 
model, we believe it to be the only tenable theory for a locally homogeneous 
surface. 

2. Line formation in an atmosphere above a backward scattering surface. 

The two models which we shall compare arc shown in Figure 1 (a) and (b). 

For a iMabert surface the scattered intensity is the sane in every direction. 

For the backscattering surface, return radiation is restricted to a narrow 
solid angle (d(0)« A real surface, with some backscattering properties, can 
be approximately represented by a superposition of the two types of surface. 

The solution for such a surface can be obtained by superposing our two solutions. 

The atmosphere in both models scatters isotropically, a close approximation 
to atomic resonant scattering. We consider only the zero phase intensity 
{ fic , iKc'^ , notation follows Chandrasekhar, I960.) 

The Intensity of radiation of a resonance line formed in a planetary atmos- 
phere is given by the solution of the radiative transfer equation. With 
standard notation, assuming plane parallel geometry, we have 
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( 1 ) 


where 

TOj . r. ti^p (- Irs^t ) 

Xc ■ Optical thickness at line center 
AXp " Doppler line width 

The incident solar beam provides an upper boundary condition: 

where fi and f*.® are cosines of zenith angles and ^ and are azimuthal 
angles; the subscript zero refers to the sun. For a lower boundary con- 
sisting of a lambert surface at » the reflected radiation is isotropic and 

the boundary condition can be expressed; 

zw f 

^ f C m' o) 

where T denotes surface reflectivity. 

For a backscattering surface, the lower boundary condition is; 

r Ij, ft,, )S; (4) 

To make the model more realistic the direct solar beam is dispersed, after 
scattering by the surface, into a narrow cone of solid angle 

The equation of radiative transfer for the two cases can be solved using 
the highly accurate approximate technique discussed In Appendix I. For TiCa) « I 
the solutions assume the following simple expressions for the emergent 


k 

intensity JL ft$) at a point on the autfacet and for the geco^trlc albedo, 
fot a spherical satellite. 


Lambert surface: 

s r/*o + 


4/t* 


( 5 ) 


f(» 




< 6 ) 


Backseat taring surface: 

( 7 ) 

fi>~) = t> - +f. - i Otr>") 

where 

fo - 

' <f u> 

For T| > 5 > I the influence of the surface is unimportant. The solutions in this 
case can be given in terms of Chandrasekhar's H*^unction 


/Fa* j' <s) 

Ptw = i C “F* 

Me can understand some features of the formation of an emission or 
absorption line from an examination of (6) and (8). Let us first discuss (6). 
Well away from a line center we observe scattered radiation, which we refer 
to as the continuum . In our models this corresponds to t, = p and the geometric 
albedo is 


PconTi 


nu um 



V 
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Close to a line center (wavelength %t the effect of a thin acattering 

layer results in the formation of a weak mission line* For a nearly 

white lambert surface, Y' I , and the contribution of the scattering layer 
to the geoaetric albedo is 

e *p ^ ^ VfC'^) 

A bright lambert surface and a layer of isotropic atomic scatterers are almost 
equally efficient (or inefficient) in scattering photons in the backward 
direction, and the atmosphere is very difficult to detect. 

The continuum geometric albedo for a backscatterlng surface is 

'f Con i ! uutn ~ T 7" 

d tv 

and, for a given value, the surface reflectively y ■ can vary from 0 to 1.0, 
depending upon the angular spread of rhe reflected sunlight, . Our so*** 
lutlon is only valid, however, fordo? small, can be large, even for a 
dark material (Y-^P), if di*? is small. Close to line center ( ^» p) the 
scattering layer contributes to the formation of an absorption line If 

a condition that will be satisfied for all T if the continuion albedo 
exceeds 0,5. The observations suggest and'^:sl : therefore 

In this case atmospheric scattering has a much greater effect on the 
emergent radiation than for a iMibcrt surface. If tl t observed continuum 
geometric albedo should be caused half by lambert scattering and half by 
backscatterlng the absorption would be nearly thirty times as effective as 
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the mission. 

In the Unit of backseat terlng alone the situation la similar to that of 
a parallel bean of light traversing a tube containing resonant scatterers* The 
parallel bean loses intensity. Other directions will gain intensity but this 
is not measured if only the parallel beam is observed. There is no limb 
darkening for the backscatterlng model and the mean air mass is two: two 
traverses of the atmosphere then lead to the factor 4 on the right hand side 
of (8). To further illustrate the differences between laanbert and backscatter- 
surfaces we show in Figure 2 (a) , (b) , and (c) a comparison of geometric 
albedo for the two cases as a function of Tj . To bring out the contrast we 
have chosen the continuum albedo at X^siC to be 0.67 In all cases. This is 
accomplished by setting V* - 1 for the lambert surface andfV* ®»67 for the 
backscatterlng surface. In the case of a backscattcring surface we have an 
additional paraneter T » proportional to the spread of the reflected 

bean. We show results for the cases T" “ 0, 0.5. For the observed values 
- 0.15, 7.0 we have Ts 0.05. Asymptotically as eo ^ ^ 

approaches the value 0.69, shown in (10). A thick atmosphere will therefore 
show line reversal. This behavior could, under sutlahle circumstances, make it 
possible to interpret an observed high resolution profile in terms of values 
for both T| and V* . 

3. A theory for a backscatterlng surface. 

We present an elementary statistical theory for multiple scattering frtmJ 
a surface consisting of microscopic reflecting facets. This is the only way 
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we have been able to identify In idilch a honogeneoua* surface can give rise to 
cero-phase gemfietric albedos in excess of unity* 

We have also considered the following possibilities: an a^rphous micro- 
structure; multiple refraction at facets; single reflection at facets. 

As»rphous micros true tore implies isotropic scattering. The geometric 
albedo for a sphere covered with an optically thick layer of conservative 
isotropic scatterecs is ^ 0.69 (F^uation 10). The question of optical 

Interference between particles (shadowing effect) is a complicated one and has 
been invoked to explain a strong phase dependence of albedo (IIapke^X963). The 
effect * is to decrease the albedo below 0.69 for non-zero phase angles 

rather than to Increase It at zero phase. 

Multiple refraction at facets can be reduced to the problem of multiple 
scattering with a phase function corresponding to the average refracted in- 
tensity for all orientations of the facets. This phase function will have <mly 
forward components, however, and a thick layer of forward scattering particles 
has a lower zero-phase geometric albedo than a layer of Isotropic scatterers. 

Finally, in Appendix 2, we offer a simple theory of scattering from a sur- 
face by single reflection from mirror facets. The mirror facets ell have the 
same zenith angle but all azimuth angles are allo^ved, consistent with the 
assumption of macroscopic homogeneity. The maxlman geometric albedo is found 
when all the mirrors point vertically. TJien we have the well-known result for 
the albedo of a polished sphere, 'p “ 0.25. 

*0ne of us had the pleasure of discussing with Professor Thomas Gold the possi- 
bility of producing the observed geometric albedo by means of surface structures, 
e.g., mountains, vertical holes, etc. We remain unconvinced that a quantitative 
theory along these lines is possible hut reserve judgement until one is presented. 
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Ue therefore consider multiple reflections not <mly as a plausible 
mechanism^ but, as far as ve can judge , the only one capable of giving zero- 
phase gecmetric albedos in excess of unity. If the directions of reflecting 
surfaces are uncorrelated the albedo is no larger th«i for sitjigle reflec- 
tions. Our model therefore Implie correlated facets or, in other \^rds, 
crystals. Fanale, Johnson, md Hatson. (1974) have .postulated the existence 
of evaporltes on the surface of lo for unrelated reas<xis. 

I^e do not assert that all scattering takes place by reflection at facets* 
The most likely source of reflectlo;ts is total internal reflection in dielectric 
crystals. Only a limited range of orientations of the crystal with respect 
to the incident beam gives rise to "cat’s eyes." Other orientations, which 
do not lead to total reflection, will exhibit scattering more similar to that 
from multiple refractions. Moreover, it is hardly likely that the surface of 
lo is uniformly covered with any single material; indeed, the reflectivity 
maps of Dollfus ^d Murray (1974) show that it is not. Me therefore picture 
a mixture of "cat's eyes" and approximations to l^bert surfaces; it is the 
theory of the former which we wish to present. 

For completeness we should draw attention to the proposal by Oetking 
(1966) that enhanced backscattering may be a diffraction effect. The lo phase 
effect shows neither the rings nor the strong polarization observed in the 
glory (van de Hulst, 1957, p. 250), and we have, therefore, not considered 
diffraction further in this paper. Our findings are, however, only applicable 
to reflecting facets much larger than the wavelength of light. 

Reflections fr<mi facets as an explanation of observed scattering properties 
has been considered by Berry (1923) , Barkas (1939) , and Middleton and 
Mungell (1952). None of the' > authors considered multiple reflections. The 
importance of multiple reflections was, to our knowledge, first pointed out 
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by Ohman (1955) but he did not present a quantitative theory. 

The facets are oriented vlth respect to the surface \^lth a statistical 
distribution Mi C/4»») » where Mw i® the cosine of the angle between the normals 
to the surfar' ■ d the facet. (See Appendix 3 for mathematical details.) The 
distribution d ^ not depend on the azimuthal angle of the facet and we exclude 
fa -ts which are oriented "Downwards,” i.e., Thus for a ray It ) to 

be k ttered into the direction in a single reflection the probability 

is proportional to where 

-Aj 

is the cosine of the normal of the facet which reflects-ff^^^-jinto^f^^^;). For 
scondary scattering we introduce a correlation function between the first facet 
and the second facet where Mr is the cosine of the angle between tlte first 

facet and the second facet. For a ray^^^^*|to be scattered into the dlrectlon{^i^^;^ 
the probability Is proportion? to 


Mw Ci,j ) s “ 
li 


2j Mt ( 

where the sum is over all possible Intermediate rays<^(/(jij^) and 


Ai'ji — 


with 


MrihUii) - « / 

^ J a C 


Ajk 


4 I 

A it)n e ) i IV^a ) Cpi ( 

Ve represent each function H| {/^») and by a gausslan function. 

. M, «xp (- (-^^J ) 

^ ) ) 

where /ij and Al^are constants to be determined by nomalization conditions. 


( 12 ) 
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The ^lergent mgles are denoted by € and ^^tbe latter refer/fa to the mtgle 
hetveen the planes containing the incident and the merr^ent rays. Each 
■catterlttg dlagr«a la nortnalised by 

.C"*** i.'**"* (U, 

In Figure 3 we sh<nf calculations ol the scattering dlagrtttts for a 
Mrrelated sirror surface for three angles of incidence and for ^ * 10* 

9^ "0*1. This ii^>lies that individual facets are distributed at random 
to the surface nonaal but that a correlation exists between pairs such that 
a related mirre * is probably to be found within - 0,1 radians. 

Hot unexpectedly, the scattering diagrams reveal a HackscatterlnR peak with 
a width of about 30*-40* superposed on a more uniform diffuse re.f lection. 

These results can be used to calculate the peoraetrlc albedo and its 
variation with solar phase, , Me assume that tin*, reflected intensity, T, 
has the simple form 

^ -f ( < 1 ^> 

where t? aid ^ are the Minnaert constants (see, for example, Veverka, 1974). 

Hotc that 5 ^ taken to be functions of cs( since this functional 

dependence is already accounted for by the scattering function for a surface 
element, ) • *^<1 a measure of limb -darkening; 

^sO and^si correspond to cases of zero and laubert lireb-darkcnlng respectively. 
6 oets the absolute reflectivity of a spherical surface, l.e., the bond 

albedo,^. For our present investigation we choose^ ■ 1 and have calculated 

* • 

two cases and^sl. TIte resulting geometric albedo and Its variation are 

calculated by the method outlined by Horrl* (1950). 'fhe gei^ctrlc albedos for 

Z 

the two cases are 1.49 and 2.00 respectively. Tlie phatte variation y " is 


II 


•ho«m 1b Figure 4, end the result indtcetes that the phase variation of lo is 
better described hy our wodel than hy a lanhert eurfaee, Wc have meager 
evidence to distinguish between the cases j^soand I since spatial resolution 
of the surface of To wuld be needed for the task. Itowever, from the 
isophotes of To taken bv Dollfus and Hurray (1974) we must conclude that lo 
does not exhibit limb-darkening » except perhaps towards the polar caps. Thus 
we would be Inclined to favor the choice of 0. 

The parameter fi'J Is critical for deterralnin?. the widtii of the bacVscattered 
peak. The width varies in proportion to ^ while the height of the peak and 
thif tero phase geometric albedo vary appicximately In inverse proportion. Thus 
ge<^etric albedos In excess of 2.H are possible, depending upon the nature 
of the surface. 

In actuality we anticipate a mixture of backscatterlnB and lamberc surfaces 
on lo, each with a reflectivity less than one. Vntll more is known about the 
surface, it Is difficult to unravel the number of parameters involved. The 
nurpose of the present calculation is to dcnonstr.ate the p/opertics of 
multiple reflections from facets, ard to shm? that this Ic n plnuslble 
model for the ohsr ed photonetric properti»»s of To. 
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Figure Captions 


Figure 1. Geometry for radiative transfer. 

In (a) the surface reflects incident radiation according to l^bert's 1^. 

Case (b) is for a backseat tering surface. 

Figure 2, Geometric albedo in the neighborhood of a line center. 

For the laidsert surface the reflectivity is chosen to be 1.0 so that the 
continuim geometric albedo is 2/3. For the backseat ter ing surface the 
continuum albedo f , = 5^ is chosen to be 2/3. We two represen- 

tative cases ( dtu) r: (0,5, 2.4); (0,0). llj is the optical thickness 
of the scattering layer at line center A=;^r (a) Z, » 3.0; (b) » 1,0; 

(c) t, = 0.1. 

Figure 3. Scattering diagrams for a surface consisting of correlated mirror 
facets. 

The angles J and are respectively the zenith angles of the incident 

and emergent rays. ^ is the azimuthal angle between the plane containing 
the Incident ray and the plane containing the emergent ray and the surface 
normal, (a), (b) and (c) illustrate the results for three angles of 
incidence: a = 6.7®, 33®, and 59®. The parameters and equal 10 and 
0,1 respectively; ^ is normalized such that J* ^ I 

Figure 4. Phase variation of the geometric albedo. 

The scattering diagram shown in Figure 3 is used for ^sO and si , aa 
explained in the text. The geometric albedos for the two cases are 1.49 
and 2.00 respectively, for conservative scattering. The result shows the 
behavior of the measured phase variation from 0® to 12® (Morrison et al., 1974), 

Figure Al. Single reflection by a. sphere c ove red with mirror facets . 

(a) Reflecting cones in the torus & i i 6 

(b) Reflection by a surface element of a cone. 

Figure A2, Schematic diagrams for multiple scattering by mirror facets. 

(a) Vector representation of incident and emergent rays, and mirror facets. 

(b) Dispersion of incident ray due to variation in ^ 

(c) Multiple scattering by mirror facets to f’th order. 


Appendix X! Calculation of emergent Int^sltv and geometric albedo. 

We solve the equation of radiative transfer (1) for IcBabert and backseat- 
terlng surfaces. The incident solar beaflB enters a scattering layer of optical 
thick <^ess *^1 at afnith angle and an asii^thal mgle (F^ure 1.). 

For simplicity vre will not carry along the index ^ in (1). Define the source 
function 




(M.) 


be expressed fomally in terms of V'C'CJ (Chandrasekhar, 
1960, p. 12) as 


^ f -rzi*ry/ct f 


X, 


^ -tS- vJ/m. 


lCx,M,i>)- I(o,M.i)e * ^ e 

At the upper boundary X^O we have 

At the lower boundary ^cTf we have, for a lambert surface, 
and for a backscatterlng surface, 

r(v„ -- j 


Jis^j (A2) 
(A3) 


(A6) 


(A5) 


(A6) 


Substituting (A2)-(A3) into (Al) and making use of (A4)-(A6) we obtain an 
integral equation for the source function 




t/ 


X 


for a Imsbert surface and 




for a backseat terlng surfacet trtiere the function Ko<x) denotes the exponential 
Integral as defined by Chandrasekhar (1960, p. 373). The Integral equations 
(A7) and (A8) can be solved iteratively on a computer. However, using a 
C(xnbination of variational and iterative techniques developed by Sze (1975), 
we can readily obtain approximate analytic solutions for the source functions 




Using (A2) we can compute the emergent Intensity In the direction ^ 

For the outer solar system we are always close to zero solar phase angle and 
the most interesting quantity is Slven by 


■Ji 

•f 


^ u. C dT f 

gir ^ ^ ( 

Jp /«• ^ 

The geometric albedo of lo can be calculated from the equation 


( '“•//‘•■’A = 


.yk, (A12) 


Appendix 2 ; Geometric albedo for a sphere covered wttli mirror facets t 
single scattering 

The only restriction comes from the assumption of local hoo»genelty. 

This requires that the mirror distribution be symmetrical with respect to 
the local vertical. Take one local vertical at a time. To maintain syim^try 
the mirrors must be formed into a cone of revolution. We shall asstme that 
the mirrors are inclined at an angle ^ with respect to the local horizon. 

If all the cones within the torus 0 i reflect backwards, the area 

of the reflecting surface, as shown in figure Al-a is 

4^ r*" ccs 0 sjn0 S$ 

The reflected light is spread . Consider the fraction of each cone 

reflecting (see figure Al-b). If all the light within the slice 2 Is 
back reflected, the fraction is he angular deviation is . 

Thus the total reflecting area is 

Ct>i $ ^ ^ ¥ 

and the total angular deviation is 

in 

The corresponding solid angle of dispersion is 

7t_ 

* f7T ( ^in O) 

From the definition of geometric albedo ve have 

^ Ct9%0 % in & 


f ( lb (if i)n 0 ') 
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It would appear that we can choose and ^ f independently. But If we 

inpose the condition that illumination is to be uniform inside the bundles we 
have 

^ i ^ 

so I 

In the limit B •s.O , we have the well-known result ^ s ^ » the geometric 

albedo for a polished sphere. 


Appendix 3: Calculation of multiple scattering by correlat; e 4 


Let us describe a mirror orientation by a unit vector m perpendicular to 
its surface, bet the incident and emergent rays be denoted by the imlt 
vectors 4j and respectively (Figure A2~a). Then the mirror which reflects 
4.J into 4|' is characterized by a vector 

aver #V 


.. 

^•j 




[^( 

If we express 4,' and 4/ in terms of angle variables and 

polar coordinates we can express /^r»(lj) » cosine of with respect 

to the vertical in the form of Ec^ . (,11) 

r m^h\) - -p- “ — 

[ .aj - fM;-/ { cPi (fl-t)) J 

For a ray that is scattered twice we have two mirrors 

-Vt' ^ J 

cosine of the angle between these two vectors is 




mil 

^ J 




where 


ic 


^ M ^ “ /■” /( „ - ( /-i/rti ^ ( l"/<n ) ^ ) 


Suppose v;e keep the incident ray fixed and vary the vector jh;j within a solid 
angle c/Wm (see Figure A2“b). It can be shovm that 


efu); = 4* 

" ^ ^ I *^j ) ) ^ 

bet the incident solar beam he In the first scattering there 

is a probability that the ray will be scattered into the upper hemisphere 

4 

as •fCM'l/fi)' ^is tay will be allowed to escape from the system ns an emergent 


2 


ray. The part that Is scattered Into the lot?er hemisphere -f will he 
scattered a second tine. The process will repeat* as shown In Fljure A2-c, «td 
the final emergent be®i the *f 

0O 

= 2*je, 

In each order of scattering we have 

IT ^ 

normalized to 
2^ f 

J 4fiH Ml C /ihi , 4 :^/ ; Mh i; ) =■ / 


The conservative nature of scattering follows from 

2tt 0 

y, 

For correlated mirror facets we introduce a correlation f unction^/^J^j, Let 


where yti7 denotes a bean in the Icwer hemisphere. Normalization is assured by 

J” /■(. ; h74i ; -«» ip ) ^1 </0‘ ip ~ ' 

Again, the system is conservative. Similarly we can compute a higher order 
scattering 


; Ml 4 p j 

=.J M Mii. ;mH i) 

To carry out the calculations numerically we reduce the radiation field to 
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a set of discrete beams. A frequently occurring Integral Is of the form 

Ue choose 20 gauss Ian points for /t and^ equally spaced points for f . The 
abo\^ Integration can be approximated by a sum 

•Si(Mh,4l)' J ^j) 

kx-i9 

\dtere we have performed a gaussian quadrature in ja and a slmpson 

quadrature in . The coefficients 4/ , hj refer to the gaussian and stepson 


weights respectively. 
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